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Abstract 

We calculate the pion self-energies in asymmetric nuclear matter in the two-loop 
approximation of chiral perturbation theory. We find three types of corrections be- 
yond the well-known linear density approximation. The resulting s-wave potential (or 
equivalent mass-shift) of a tt~ in the center of a heavy nucleus like Pb turns out to 
be = Am^- ~ 14 MeV, about half of what is needed to form the deeply bound 
and narrow pionic atom-states recently observed at GSI. The potential for a vr"*" under 
the same conditions is = Am^+ ~ — IMeV. In the same way we calculate the 
mass-shifts of "''-mesons in symmetric nuclear matter and neutron matter and find 
an increase of the K^-mass by 9% and 5%, respectively. As a further application we 
give an analytical result for the third order term in the scattering length expansion of 
the in-medium self-energy of an interacting light-boson/heavy- fermion system. 
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The modification of pion (or kaon) properties in dense nuclear matter is characterized 
by their density and four-momentum dependent self-energy. Systematic studies of pion- 
nucleus scattering and pionic atom data have lead to a standard parametrization of the 
corresponding phenomenological pion-nucleus optical potential (for a review see [ij]). On 
that basis Friedman and Soff and Toki and Yamazaki 0, |^ predicted the existence of 
deeply bound and surprisingly narrow pionic atom-states in heavy nuclei. Pionic Is and 2p 
states in ^°^Pb and ^osp]-, i^ave recently been observed at GSI p, ^ using (d,^He) transfer 
reactions, and there is on-going experimental activity to investigate more detailed aspects of 
such states like systematic isotope-shift effects. These deeply bound pionic states owe their 
existence to a subtle balance between the attractive Coulomb potential and the repulsive 
s-wave optical potential. The net attraction is localized at and beyond the nuclear surface. 
The resulting pion wave function has little overlap with the nucleus and this explains the 
observed small absorption widths, P < 1 MeV. Our aim here is to provide a reliable baseline, 
guided by chiral symmetry, for the discussion of these phenomena. 
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The starting point in the construction of the pion-nucleus optical potential is the low- 
density theorem [0, |^. It states that, to linear order in the proton and neutron densities, the 
optical potential is given by the respective TriV-scattering threshold T-matrices times these 
densities (see eq.(4) below). Corrections beyond the linear density approximation have been 
derived by Ericson and Ericson ^ for symmetric nuclear matter employing analogies with 
optics. The correction term of can be interpreted as double-scattering in the nuclear 
medium, with the relevant vrA^-interaction parametrized through the scattering lengths. 

Deeply bound pionic states are particularly sensitive to the s-wave part of the pion-nucleus 
optical potential P]. The s-wave vrA^-interaction is, on the other hand, strongly constrained 
by chiral low-energy theorems (i.e. the Tomozawa- Weinberg relation for the vrA^-scattering 
lengths). The tool to systematically investigate the consequences of spontaneous and explicit 
chiral symmetry breaking in QCD is chiral perturbation theory. Observables are calculated 
with the help of an effective field theory formulated in terms of the Goldstone bosons (vr, K, rf) 
and the low-lying baryons. A systematic expansion in small external momenta and meson 
masses is possible. It is the purpose of this work to calculate the s-wave pion (and kaon) self- 
energies in asymmetric nuclear matter completely up to sixth order in the small momentum 
expansion. The results include (to that order) all corrections beyond the linear density 
approximation generated by the chiral meson-baryon interactions together with their explicit 
dependence on the proton and neutron densities. 
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Fig.l: Some in-medium self-energy diagrams up to two loops. Solid and dashed lines repre- 
sent nucleons and pions, respectively. The combinatoric factor of the last diagram is 1/2. 



From the start we will restrict the calculation of the in-medium self-energies to on- 
shell pions (and kaons) with zero three-momentum, i.e. = {rriT^K,^)- The only new 
ingredient to perform calculations at finite baryon density (instead of scattering processes 
in the vacuum) is the particle-hole or in-medium nucleon propagator. It expresses the fact 
that the ground-state has changed from an empty vacuum to a filled Fermi-sea of nucleons. 
In the heavy baryon limit the in-medium nucleon propagator reads 



Po + 
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with the energy pq counted modulo the large nucleon mass M and kp^n are the Fermi- 



momenta of the protons and neutrons, related to their densities Pp^n = ^p,n/37r^ in the usual 
way. With the help of the projection operators (1 ± T3)/2 one can easily treat the case 
of different proton and neutron densities (i.e. isospin-asymmetric nuclear matter). Note 
that the in-medium propagator eq.(l) splits additively into the vacuum propagator and a 
medium-insertion. This allows to organize the diagrammatic calculation according to the 
number of medium insertions. Diagrams with no medium insertion just renormalize the 
meson mass to its (measured) vacuum value. Diagrams with one medium insertion (such as 



the first one in fig. 1) produce exactly one factor of density Pp or p„ and thus contribute to 
the hnear density approximation. Corrections beyond the hnear density approximation arise 
only from diagrams with (at least) two medium insertions (see fig. 1). Their evaluation leads 
to an integral over two Fermi-spheres of different radii which can be reduced to a simple 
one-parameter integral by the following formula, 
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with the weight-function 
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Furthermore, working at g* = simplifies the calculation considerably since all two-loop 
diagrams in which the in- or outgoing pion couples directly to the nucleon vanish identically 
{a ■ q = 0). The third diagram in fig. 1, with three nucleon propagators, requires some 
special attention. Here, one encounters the square of the in-medium propagator eq.(l) 
whose medium-insertion turns out to be proportional to the derivative of a delta-function 
S'{po). This together with the fact that the internal pion-propagator is an even function of 
Po, implies that there is no correction beyond the linear density approximation from this 
particular diagram. 

We are now in the position to present the complete two-loop result for the pion self- 
energies in asymmetric nuclear matter. In the case of a negatively charged vr" it reads: 



(/Cp, kn) 



kl 



37r2 



kl r 



37r2 



+ n 



relikpi kn) 



+ Il^orikp,kn). (4) 



The first two terms represent the linear density approximation with = AiT{l+mT^/M)a^ , 
the isospin-odd and isosp in-even vrA^-scattering threshold T-matrices calculated to 1-loop 
order in chiral perturbation theory |]TU|, |TT|, 
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Here, /tt = 92.4 MeV is the weak pion decay constant and qa = 1-27 the nucleon axial- vector 
coupling constant. The combinations of low-energy constants C = C2 + — 2ci and 6(A) 



have been discussed in detail in refs. [p!0| , |TT[ and the vanishing of the m^-term for T~j^ has 
been proven in ref . |Tl| . Experimental values for have been extracted from the shift and 
width of the Is level in pionic hydrogen by the PSI group |]T2[, with the final result 



1.847 ± 0.086) fm, 
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It is important to note that the pion-loop correction closes the gap between the Tomozawa- 
Weinberg prediction 771^^/2 = 1.61 fm and the empirical value of T~^. In the case of 
TJ^ no accurate chiral prediction is possible since the next-to-leading order pion loop term 
3g\m^/6ATTf^ = 0.18 fm is still large in comparison to the (anomalously) small empirical 
value of TJ]v (which is actually compatible with zero) . 



The last two terms in eq.(4) comprise the corrections beyond the hnear density approxi- 
mation. The first term 

is a small relativistic correction from the second diagram in fig. 1. It originates from the 
energy-dependent pseudovector nN Born amplitude in forward direction, integrated over 
the Fermi-sphere. More interesting are the contributions from the last three diagrams in 
fig. 1, with two medium insertions. They describe scattering processes in which two nucleons 
in the Fermi-sea participate, with the result 
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The first line in eq.(9) is the Ericson-Ericson double scattering term generalized to asym- 
metric nuclear matter {kp 7^ A;„). The term proportional to g\ is new. It represents the 
effect induced by vrTr-interaction with two virtual pions being absorbed on the nucleons in 
the Fermi-sea. Note that the off-shell 47r-vertex is not uniquely defined by the chiral La- 
grangian and therefore one should view the last two diagrams in fig. 1 always as one entity. 
Such ambiguities hint at the fact that an off-shell in-medium self-energy itself is not a unique 
concept in field theory. Interestingly, the two-nucleon correlation term n~^(fcp, /c„) is an even 
function of the proton-neutron asymmetry. 

Next, we turn to the vr"*" self-energy in asymmetric nuclear matter. It can be easily 
obtained from the vr" one by an isospin transformation which exchanges the role of protons 
and neutrons, 

TV^{kp,kn) = Tl-{kn,kp). (10) 
For the sake of completeness we also give the result for the (neutral) vr'' in-medium self-energy, 

n°(A;p, kn) = —T^N ^2^2 ■'^cor(^P' ^n) , (H) 

with the corresponding two-nucleon correlation term. 
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The differences in comparison with the charged pion case can be traced baclc to certain 
different isospin factors. 

The self-energies are related to the threshold s-wave optical potentials C/^°(/Cp, by 
IP^{kp, kn) = 2771^^11'^^ {kp, kn), and we can identify in-medium pion mass-shifts, Am^^o = 
U^'^{kp, kn) in order to give a physical interpretation to these potentials. For a numerical 
example we consider asymmetric nuclear matter of density {kp + /c^)/37r^ = 0.165 fm~^ with 
a neutron-to-proton ratio N/Z = {kn/kpY = 1.5 as it is typical for a heavy nucleus like Pb. 
This implies proton and neutron Fermi-momenta of kp = 246.7 MeV and fc„ = 282.4 MeV. 
The in-medium pion mass-shifts Arrij^^o are found in such a dense environment as 

Am^- = (9.5 - 1.0 + 6.0 - 0.7) MeV = 13.8 MeV , (13) 

Am^+ = (-7.5 + 1.0 + 6.0 - 0.7) MeV = -1.2 MeV , (14) 

Am^o = (1.0 + 0.0 + 5.7 - 0.6) MeV = 6.1 MeV . (15) 

Here, the first number in brackets gives the contribution from the linear density approxi- 
mation using the central PSI- values for (see eq.(7)). The second number corresponds 
to the small relativistic correction eq.(8) (zero for 7r°), the third number gives the double 
scattering contribution and the fourth number corresponds to the new g\-teTms induced by 
TTTT-interaction. Note that the splitting of vr" and 7r+ in asymmetric nuclear matter comes al- 
most entirely from the leading order T^-term driven by the Tomozawa- Weinberg low-energy 
theorem. This splitting is absent in symmetric nuclear matter. Clearly the Ericson-Ericson 
double scattering term ^ is the dominant correction beyond the low density approximation 
and it is almost equal for charged and neutral pions even in asymmetric nuclear matter. The 
effect induced by yrvr-scattering reduced this by about 10%. Both these two-nucleon correla- 
tion effects show up in a calculation of the pion-deuteron scattering length There, the 
^f^-term contributes only at the 3 — 4% level due to the much smaller densities relevant in 
the deuteron. 

A repulsive mass-shift Am^r- — 14 MeV is still too small (by about a factor of two) to 
form the deeply bound pionic atom-states in Pb as narrow as they are observed. We have 
therefore calculated some contributions of seventh order in the small momenta m.„, kp, k^. 
For example, continuing the 1/M-expansion of the second diagram in fig. 1 one obtains 

SgWM + f'n) .... 

n^^l{Kp, Kn) - 407r2M3 /2 ' ^ ' 

which results in an additional small downward mass-shift of Am^r- = —0.3 MeV. In the 
dominant double scattering diagram (fourth graph in fig. 1) one can either correct for the 
nucleon kinetic energies in the pion-propagator or replace the upper or lower leading order 
vertex by a second order vertex proportional to C — g\/8M. Both types of corrections give 
together, 

n-.(^p, A;„) = ^^§§^{Kkn + likl + kl)[\n + g\- 8MC] } , (17) 

and numerically they cancel each other, Am^- = (—0.2 -|- 0.2) MeV. 



In phenomenological calculations [|1^] it has been common practice to introduce a term, 
— 167r_Bo PpPn, in the tt" self-energy and to choose Re So such that the missing repulsion is 
accounted for. A term of this form would show up at three-loop order in connection with 
pion-absorption, nNN NN. This process involves large momenta of order ^fMml which 
are difficult to cope with in the counting scheme of chiral perturbation theory (see also the 



discussion of the inverse process A^A^ — > NNn at threshold in ref. |T^). A large negative 
value of Rei?o [Q would be at odds with theoretical many-body calculations [16], though 
within considerable uncertainties [O. An alternative way to generate the missing s-wave 



repulsion could be the in-medium renormalization of the pion decay constant f^^ 

The two-loop ChPT calculation of the pion self-energy in asymmetric nuclear matter can 
be easily taken over to the case of kaons [K^, K^). The mass-shifts of anti-kaons {K~ , K^), 
on the other hand, are strongly affected by the (non-perturbative) in-medium dynamics of the 
subthreshold A(1405)-resonance At two-loop order the self-energy in asymmetric 
nuclear matter reads 
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matrices. The empirical values T^\; = 0.4 fm and T^j^j = —6.3 fm stem from a dispersion 
relation analysis of KN and ^A^-scattering by Martin |^ . The expression in the second line 
of eq.(19) derives from the double scattering graph (fourth diagram in fig. 1) using chiral KN- 
interaction vertices. Here we have introduced the weak kaon decay constant /i^- = 113 MeV in 
order to have a more realistic leading order value for T^j^ = —frtK/ fx = —7.6 fm. Of course, 
we have also evaluated the corrections beyond the linear density approximation which come 
from the 1/M-expansion of the pseudovector KN Born-amplitude and from the diagrams 
with virtual (vrvr , tt^t] , T^r^) -exchange between the kaon and two nucleons in the Fermi-sea. 
We found that these lead to negligibly small kaon mass-shift of 0.6 MeV or less. Therefore 
we have omitted the corresponding expressions in eq.(19). 

The K^ mass-shift in symmetric nuclear matter of density 0.165 fm~'^ is ArriK = (30 + 
15) MeV = 45 MeV according to eq.(19), with the first and second number coming from the 
linear density approximation and the double scattering term, respectively. In pure neutron 
matter at the same density the K^ mass-shift is substantially smaller, ArriK = (19 + 6) MeV 
= 25 MeV. The coupled channel calculation of |]19[ has found very similar results. Again, the 
K^ self-energy can be immediately obtained by an isospin transformation as IlKo{kp, k^) = 
^K+{kn,kp). In this context, we have to comment on the "model independent" result for 



.(1) 



the K~^ mass-shift of Lutz [^. For symmetric matter {kp = k^) (and to the order in 
itik/M considered here) we agree with the real part of his double scattering term. The 
claimed imaginary part (implying a K^ decay width of Tk+ = 4 MeV) is however wrong. 
For a K^ at rest all on-shell scattering processes to change any of its quantum numbers are 
Pauli-blocked. 



b 6 o o 



Fig. 2: In-medium self-energy diagrams up to third order in the scattering length a. 



As a side remark and a further application of effective field theory in the same context, 
we consider the self-energy a light scalar boson moving in a medium consisting of heavy 
spin-1/2 fermions. The boson-fermion interaction is assumed to be completely described by 
the scattering length a, i.e. a local contact interaction of the form £j„j = 2TTaHH(fP'. From 
the diagrams in fig. 2 one can compute the in-medium self-energy n(a, kf) up to third order 
in the scattering length a, with the analytical result 

n(a, kf) = -^ak) + -^a'k)- ^-^f^ a' kj , p = ^ . (19) 

The new triple scattering term (which goes beyond the work of ref.0) originates from the 
two topologically distinct three-loop diagrams in fig. 2, either with three medium insertions 
or with two medium insertions at the lower and upper fermion-propagator. The latter 
combination is not a renormalization of the a'^kj double scattering term. Eq.(20) cannot be 
directly applied to the irN and KN systems due to the additional isospin quantum number 
and the host of momentum dependent chiral meson-baryon vertices. Nevertheless, the result 
eq.(20) may be useful for certain condensed matter systems. 

In summary, we have used chiral perturbation theory to calculate the complete pion and 
kaon self-energies in asymmetric nuclear matter up to two-loop order. We have found new 
corrections beyond the linear density approximation together with their explicit dependence 
on the proton and neutron Fermi-momenta. To that order the n~ mass-shift in Pb is about 
half of what is needed to form the observed deeply bound and narrow pionic states. The 
present result can also be of some guidance for the analysis of vr"*" electroproduction off ^He 
P^ . For the mass-shift we confirm earlier results of model calculations. 
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